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Abstract 

The  spherical  lattice  gas  is  a  modification  of  the  ordinary  lattice  gas 
in  which  the  occupation  mamber  of  each  cell  is  permitted  to  be  any  real 
number  rather  than  +  1.   However  the  sum  of  squares  of  the  occupation  numbers 
is  required  to  equal  the  number  of  cells.  This  permits  one  to  evaluate  the 
partition  function  by  integrating  over  the  surface  of  a  certain  sphere  rather 
than  by  summing  over  lattice  points  on  that  surface.  The  partition  function  and 
the  equation  of  state  of  the  gas  are  evaluated  in  this  way.   It  is  foimd  that  in 
three  dimensions  the  gas  condenses,  but  not  in  one  or  two  dimensions.   Graphs  of 
the  phase  transition  curve  and  of  the  isotherms  in  three,  two  and  one  dimensions 
are  presented. 

The  analytical  work  is  simplified  by  taking  advantage  of  the  relationship 
between  the  properties  of  the  lattice  gas  and  of  the  Ising  model  of  a  ferromagnet. 
This  relationship,  demonstrated  by  C.N.  Yang  and  T.D.  Lee'-  -'  for  the  ordinary 
lattice  gas  and  Ising  model,  also  applies  to  the  spherical  lattice  gas  and  the 

spherical  model  of  a  ferromagnet.   The  properties  of  the  latter  have  been  eval- 

\2\ 
uated  by  T.H.  Berlin  and  M.  Kac ^  ^ .   Graphs  of  the  Isotherms  of  the  spherical 

model  of  the  magnet,  which  were  found  in  the  course  of  the  work,  are  also  presented. 
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1.   Introduction. 

The  spherical  lattice  gas  is  a  simplified  model  of  a  gas.   We  introduce 
it  because  we  can  deduce  its  equation  of  state  statistical  mechanically.   We 
find  that  it  undergoes  phase  transition  in  three  dimensions,  although  not  in 
one  or  two  dimensions.   We  also  find  that  its  isotherms  and  transition  curve, 
shown  in  Figure  9  are  in  qualitative  agreement  with  experimental  curves . 
However  they  exhibit  a  physically  unreasonable  negative  pressure  at  very  large 
and  very  small  specific  volumes.   This  can  be  traced  to  an  inadequacy  of  the 
model  for  such  specific  volumes. 

The  ordinary  lattice  gas  is  defined  by  first  dividing  the  volvmie  V,  which 
the  gas  occupies,  into  V  identical  cells  (e.g.  cubes),  each  of  unit  volume. 
The  potential  energy  of  interaction  between  two  molecules  of  the  gas  is  defined 
to  be  infinite  if  they  occupy  the  same  cell,  to  be  -kj   <  0  if  they  occupy  cells 
which  are  nearest  neighbors,  and  to  be  zero  otherwise.   Thus  at  most  one  particle 
can  occupy  each  cell.   Therefore  the  distribution  of  particles  can  be  described 
by  the  V  coordinates  c  ,  CT  ,  ...,cr  where  cf.    =   +1  if  cell  i  is  occupied  and 
6.    =   -1  if  it  is  not.   We  define  the  spherical  lattice  gas  by  permitting  each 
cr.  to  have  any  real  value,  but  requiring  that  the  a.    lie  on  the  sphere 

r  ^i  =  V  (1) 

1=1 

Then  the  grand  partition  function  Q   of  the  spherical  lattice  gas  is  obtained 
from  that  for  the  ordinary  lattice  gas  by  replacing  sums  over  lattice  points 
which  lie  on  the  sphere  (l)  by  integrals  over  that  sphere. 

We  were  led  to  Introduce  the  spherical  lattice  gas  by  the  work  of  C.  W.  Yang 
and  T.  D.  Lee  ^ -I  and  that  of  T.  H.  Berlin  and  M.  Kac^J.   Yang  and  Lee  showed 
that  Q  ,  the  grand  partition  function  of  the  ordinary  lattice  gas,  could  be 
expressed  in  terms  of  Q  ,  the  partition  function  of  the   Islng  model  of 


dimensional  magnet  nor  for  a  two  dimensional  one  with  H  j^  0^  although  the  spontan- 
eous magnetization  is  known  as  a  function  of  temperature  in  the  two  dimensional 
case.   Therefore  the  phase  transition  curve  of  the  two  dimensional  lattice  gas 
could  be  determined,  but  not  the  isotherms,  while  nothing  could  be  found  for  the 
three  dimensional  case. 

Berlin  and  Kac  introduced  the  spherical  Ising  model,  which  is  related  to  the 
ordinary  Ising  model  as  the  spherical  lattice  gas  is  to  the  ordinary  lattice  gas. 
They  were  able  to  evaluate  its  partition  function  Q   for  any  value  of  H  in  one, 
two  and  three  dimensions.   It  exhibits  spontaneous  magnetization  in  three  dimensions 
but  not  in  one  or  two  dimensions. 

We  noticed  that  the  method  of  Berlin  and  Kac  could  be  used  to  evaluate  the 

grand  partition  function  of  the  spherical  lattice  gas.   We  also  noticed  that 

Q^^  is  related  to  Q„^  in  the  same  way  that  Q^  is  related  to  Q^.   Therefore  we 
dG  ol  (j  1 

could  actually  use  their  results  to  obtain  the  value  of  Q^.,. 

In  the  course  of  our  analysis  of  the  isotherms  and  transition  curve  of 
the  spherical  lattice  gas,  we  were  led  to  examine  the  magnetization  curves  and 
spontaneous  magnetization  of  the  spherical  Ising  model.   Therefore  these  curves 
are  also  included. 

2.   Basic  formulas. 

The  grand  partition  function  of  the  lattice  gas  is  [3,p.384J 

QG(V,;,y)  =  i:  ^^••-  11  ^^  exp[|  YZ_^^   a^.^^crj  exp[(cy.  |  In  y)      (2) 


V 
X  1 

j 


-aj  exp[i(c^.lny)v] 


Here  6.   equals  +1  if  cell  i  is  occupied  by  a  particle  and  -1  if  it  is  not. 
We  note  that  the  a.    satisfy  (l)  and 

af  =  1  .  (5) 

The  constant  fi-  is  defined  in  terms  of  Boltzraann's  constant  k  and  the  absolute 
temperature  T  by 

^    =  j/kT  .  W 

The  quantity  a. .  equals  +1  if  cells  i  and  j  are  nearest  neighbors  and  is  zero 
otherwise,  while  c  is  the  number  of  nearest  neighbors  of  a  cell.   For  cubic  cells 
in  one,  two,  and  three  dimensions  c  equals  two,  foui)  and  six  respectively.   The 
fugacity  is  denoted  by  y. 

From  the  definition  of  the  spherical  lattice  gas,  its  grand  partition 
function  is 

V 


Qg^(V,^,y)=       j---^     da^...da-^exp[i-|Z_^a..a.aJexp[(c^.|lny)         (5) 


X 


y~     ^2    expft   (cy+  In  y)vj     . 
The  partition  function  for  the   Is-ing  model  of  ferromagnetism  is     [3,p.353] 

Q.(nJ,?/)    =X:     ---II        expfZZ       ^^fi<^.^   ^      Her]     .  (6) 

In  (6)  N  is  the  total  number  of  lattice  sites  and 

H  H 

^'w     ■  '^' 

Here  ^   is  the  magnetic  moment  per  particle. 


-  k 


The  partition  function  for  the  spherical  model  of  ferromagnetism  is  [2;, p. 827] 

(I  rQ   ^  N   ^ 

Q3,(»,^,W)=  J  ...|  dcr^...dcv„exp[f  E^ay«.«.  .>/Z  <.jj  .      (8) 

E  Oi  - » 

In  (2)  let  us  replace  V  by  N  and  set 

In  y  =  2[^-  c^]  . 
Then  by  comparing  (2)  with  (6)  we  obtain  the  following  relationship^,  first 
obtained  by  Yang  , 

^^U(^,e^v[2{n-   ^J^)])=  exp[?f  -  I  c^J  Q^(nJ,M)  .  (9) 

Upon  comparing  (5)  and  (8)  we  obtain  the  corresponding  relation  for  the 
spherical  models, 

Qg/N,(^,exp[2(  H-    cj^)]^  =  exp[;^  "  I  =  ^J  %^i^,(^,n)     •  (lO) 

From  either  (9)  or  (lO)  it  follows  that  [j^p.jQ^J  the  thermodynamic  properties 
of  the  lattice  gas  are  related  to  those  of  the  magnet.   These  relations  are  the 
same  for  the  spherical  and  the  original  models.   They  are 

p*((>,y)  =  ^-^  - 1  ,  (11) 

A},V)-^^^     .  (12) 

Here  p  =  p/J  where  p  is  the  pressure  of  the  gas,  v  is  its  specific  volume,  and 
%  =   M/i-i  where  M  is  the  magnetization  of  the  magnet.  F,  the  limiting  free  energy 
per  particle  of  the  magnet,  is  defined  by 

N->oo 


Here  Q  denotes  either  Q  or  Q^y-  TH    is  given  by 

Equations  (ll)  and  (l2)  are  the  parametric  form  of  the  equation  of  state  of 
either  the  lattice  gas  or  the  spherical  lattice  gas.   The  parameter  is  fir, 
which  ranges  from  minus  infinity  to  plus  infinity. 

To  obtain  the  equation  of  state  of  the  spherical  lattice  gas  we  must  insert 
^{'^  >  ^  ))    "the  limiting  free  energy  per  particle  of  the  spherical  Ising  model,  into 

(ll)  and  (l4)  and  then  use  (l^)  in  (12).   For  a  simple  cubic  lattice  in  n  dimensions 

["2] 
(n  =  1,2,3),  Berlin  and  Kac ^  J  have  evaluated  this  free  energy  which  we  denote  by 

F^iV  ,f   )  .   They  obtain 

F„(W,{?)=iln(|:)-if„(z).^..ifJii^.  (15) 

In  (15)  we  have  omitted  an  additive  constant  which  occurs  in  [2]  due  to 
multiplication  of  Q   by  a  normalization  factor.   Here  z  is  a  real  positive 
solution  of  the  equation 


^2        dfjz) 


2l3(z-n)2 


(16) 


The  function  f  (z)  is  defined  by 


2rt   2ir 


^J^^='^        f.  ..   dw    dw  mfz  -  X:  cos  w  I  .  (17) 

It  will  prove  convenient  later  to  note  that  (16)  implies  dF  (z)/dz  =  0  and  that 
z  is  an  even  function  of  "H .      Then  (I5)  shows  that  F   is  also  an  even  function  of 

u. 

The  equation  of  state  of  the  spherical  Ising  model  of  a  magnet  is  obtained 
by  inserting  (I5)  into  (l^).   It  is 


1.00 


7n  .50 


Fig.  1  -  Isotherms  of  the  one-dimensional  Ising  model  of  a  f erromagnet .   The  curves, 
based  on  equation  (20),  give  the  normalized  magnetization  Jn   =  M/m-^  as  a 
function  of  the  normalized  external  field  ^   =   H/kT  for  constant  values  of 
(J  =  j/kT.   Since  ^   is  an  odd  function  of  -V  ,  only  the  non-negative  values 
of  A  are  shown. 


From  (18)  we  see  that  %   is  an  odd  function  of  ?^  .   In  the  following  sections 
we  shall  examine  the  equation  of  state  of  the  spherical  lattice  gas  and  of  the 
spherical  Ising  model  in  detail  on  the  basis  of  the  above  equations. 

In  order  to  compare  the  properties  of  the  spherical  lattice  gas  with  those 
of  the  lattice  gas,  we  must  determine  the  properties  of  the  latter.   We  can  do 
this  in  the  one  dimensional  case.  The  limiting  free  energy  per  particle  of  the 
one  dimensional  Ising  model  F  ,  is  [3, p. 585] 


F^jCN,^  )  =  m  e^  cosh  V+   (e2^sinh^?^+  e'^l) 


1/2 


(19) 


Upon  substituting  (19)  into  (l^)  we  obtain 


%. 


1/2 


^^{'H,  f)    =      (sinh^y-    +  e"^<?)     sinh  ?^ 


(20) 


(This  formula  is  misprinted  in  [3]-)   From  (20)  we  deduce,  omitting  the  subscripts 

'niiV-,0)    =  tanh  'H 

7/l(W-,co)  =1   7^  >  0 

7ll{0,^   )  =  0    y  <  CO 

^(co,^  )  =  1  . 
Equation  (21c)  demonstrates  the  well  known  fact  that  spontaneous  magnetization 
does  not  occur  in  the  one-dimensional  Ising  model.  Magnetization  curves  based 
on  (20)  are  shown  in  Figure  1.   Less  complete  graphs  are  given  in  [3]. 

Upon  substituting  (19)  into  (ll)  and  (20)  into  (12)  we  obtain  the  equation 
of  state  of  the  one  dimensional  lattice  gas  in  parametric  form 

1/2 


^  imp  cosh  -^  .   (/^  sinh^  ^  .   .-^^y  .    A 


v   =  2[l.      (s 


sinh   y    +   e      /  )  sinh  n 


(21a) 
(21b) 
(21c) 
(21d) 


(22) 


(23) 


6a  - 
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Isotherms  of  the  one  -  dimensional  lattice  gas  based  on  (22)  and  (23)  are  shovm 
in  Figure  2.   We  see  from  (23)  that  as  yj-   varies  from  +  co  to  -  oo,  v  varies 
continuously  and  hence  no  phase  transition  occurs.   We  also  obtain  easily  the 
following  values,  listed  in  Table  I 

Table  1 

"H-  7?L  V  P  ■ 

2  H 

+  (X)  +1  1  -^  —7 >  +  CO 


/ 


0  0 

-  oo  -  1  +00 


2  -^   ln(2  cosh  f   )  - 


We  see  that  the  equations  of  state  of  the  one-dimensional  Ising  model  and  of 
the  one  dimensional  lattice  gas  are  physically  quite  reasonable. 

3.   The  one- dimensional  spherical  lattice  gas. 

Since  the  properties  of  the  spherical  lattice  gas  can  be  found  in  terms  of 
those  of  the  spherical  Ising  model,  we  shall  first  examine  that  model.   Its 
equation  of  state  is  given  by  (l6)  and  (l8)  which  become,  in  one  dimension, 

'^    2(iiz-ir     (z -1)  -^ 


(24) 


From  {2k)   and  (25)  we  find  that 

772(7^,0)    =  2?/[l   +   (l+i.;V2)l/2j      ,  (26a) 

%iU  ,oo)    =   1   ,                          y  >   0   ,  (26b) 

7n(o>  Q  )  =  0    ,                   i^  <  o°^  (26c) 

Uioo,}.)    =  1  ■      (26d) 


-  7a 


1.00 


7k   .50 


Fig.  3"  Isotherms  of  the  one-dimensional  spherical  model  of  a  ferromagnet  based  on 
equations  (2k)   and  (25).   The  normalized  magnetization  97P  =  M/n  is  shown 
as  a  function  of  the  normalized  external  field  y    =  H/kT  for  various  con- 
stant values  of  y  =  j/kT.   These  curves  lie  below  the  corresponding  ones 
of  the  Ising  model  shown  in  Figure  1. 


-  8  - 


The  proofs  of  these  results  are  similar  to,  but  simpler  than,  those  given  in 
Appendix  B  for  the  three  dimensional  case.   Equation  (26c)  shows  that  spontaneous 
magnetization  does  not  occur  in  this  case.   Magnetization  curves  of  the  spherical 
Ising  model,  based  on  {2k)   and  (25),  are  shovm  in  Figure  5-   For  a  given  fj-   and  (j   , 
the  magnetization  of  the  spherical  model  is  less  than  that  of  the  Ising  model 
(Figure  l) . 

From  (11)  and  (15)  we  obtain  for  the  spherical  lattice  gas  in  one  dimension 

F  +  ^ 
i{(^,n)    =^ -1   ,  (28) 

""z  .  (z^-l)l/2 


F^(>^,  j^)  =|ln(j)  -  |ln| 


^  (29) 


4a(z-i) 


From  (12)  and  (25)  we  obtain 

(30) 


^^}'^^-'^^2fi^ir 


Equations  (28)  -  (jO)  express  the  equation  of  state  of  the  one  dimensional 
spherical  lattice  gas  in  temis  of  the  parameter  ?r  .      The  quantity  z  is  defined  as 
a  root  of  (2^).   Isotherms  of  the  one-dimensional  spherical  lattice  gas  based  on 
(28)  -  (50)  are  shown  in  Figure  h.      The  method  of  computation  used  in  preparing 
this  and  the  subsequent  figures  is  explained  in  Appendix  A.   Since  spontaneous 
magnetization  does  not  occur  in  the  spherical  magnet,  phase  transition  does  not 
occur  in  the  spherical  lattice  gas.   Much  as  in  Appendix  B  we  obtain  the  special 
results  shown  in  Table  2  • 


3       CD      ^^ 


Table  2 


7)v 


2%        ^ 
^   y   ->  +  00 


"    2  "/2" 
1+1+ fl^) 


We  note  from  the  third  line  of  the  table  that  when  v  becomes  large,  p  becomes 
negative.  Furthermore  for  large  v  the  isotherms  corresponding  to  large  values 
of  y  lie  below  those  corresponding  to  smaller  values  of  j?  .   This  unphysical 
behavior  of  the  equation  of  state  for  large  v  is  due  to  the  inadequacy  of  the 
spherical  model  in  this  limit. 

h.      The  two-dimensional  spherical  lattice  gas. 

For  n  =  2,  equations (l6)  and  (l8)  for  the  spherical  Islng  model  become 

^^   -      /'  2  ^  f^^(f)'  (51) 

'    2a{z-2f  ""^        ^ 

In  (31),  K(a)  is  the  complete  elliptic  integral  of  the  first  kind,  0  <  a  <  1. 
From  (31)  and  (32)  we  find  that  (26)  is  still  valid  for  n  =  2.   This  is  shown 
as  in  Appendix  B  by  using  the  facts  that  K(0)  =  Jt/2  and  K(l)  =  oo .   From  (26c) 
we  see  that  spontaneous  magnetization  does  not  occur  for  the  spherical  Ising 
model  when  n  =  2.   The  Ising  model,  however,  does  exhibit  spontaneous  magnetization 
when  n  =  2.   Magnetization  curves  based  on  (31)  and  (32)  are  shown  in  Figure  5. 


1.00 


"K  .50 


3.0 


Fig. 5  -  Isotherms  of  the  two-dimensional  spherical  model  of  a  ferromagnet  "based 
on  equations  (jl)  and  (52).   The  normalized  magnetization  %,  -   m/ii  is 
shown  as  a  function  of  the  normalized  external  field  ,4-'  =  H/kT  for  various 
constant  values  of  (/  =  j/kT. 
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To  obtain  the  equation  of  state  of  the   spherical  lattice   gas   in  two 
Limensions,   we   set  n  =  2   in   (ll)   and   (15)   which  become 
F^  +  % 


Here  fp(z)    is  defined  by 


^2^^)    =  2  dw^dw^  Infz   -   cos  w^   -   cos  wH    . 

^    "^^        00 

Upon  integrating  once,  fp(z)  becomes 

f2(z)  =\        j    dw  In  I  j(z  -  cos  w)  +  [(z  -  cos  w)^  -  l]    >    .  (36) 

o 

From  (12)  ani  (52)  we  obtain 

Isotherms  of  the  two-dimensional  spherical  lattice  gas  based  on  (3^)-(57) 
are  shown  in  Figiire  6.  From  (26),  (3^)  and  {jj),   by  considerations  like  those 
in  Appendix  B,  we  obtain  the  special  results  in  Table  3. 


Table  3 


^ 


P 


2^ 


f  |-^2^^m)^y4 


In  1^1 
-1  +00  -     \       I    -> 


«^  -o 


-  11  - 

Here  z  is  the  solution  of  equation  (5)  in  Appendix  C.   We  see  from  the  third 
line  of  the  table  that  for  large  v  the  behavior  of  the  two  dimensional  spherical 
lattice  gas  is  unphysical,  just  as  in  the  one  dimensional  case. 


5.   The  three-dimensional  spherical  lattice  gas. 

In  three  dimensions  (16)  and  (18)  for  the  spherical  Ising  model  become 


2«(2-3)' 


(38) 


^^V,p-jp-)     ■ 


In   (38)    df   /dz   is   given  by 


df^Cz) 


2n  2rt        27t 


dz 


dw  dw  dw 


/^   w        ]      I  /      z   -   cos  w     -   cos  w     -   cos  w 

^2^^  000  ^  ^  ^ 


(39) 


(^0) 


^       (       du(„     ^_  J  K(„   _   L,  ,J 
If  ' 


Z     -     cos    W  'z     -     COS    W 


In  Appendix  B  we  show  that 


U(7^>o)  =  2?/  1  +  (1  +  k)^^) 


1/2 


n-l 


(41a) 


n,oo)   =  1 


1/2 


(0+,/)    =    {1   -  ^) 


[0, 


p-0 


n>  0 

for     0    >   i=    .2527. 

;<;c 


(iiib) 

(iHc) 

(4ld) 


%{o,,(j)    =  1  .  (i^le) 

Equations  (4lc)  and  (4ld)  show  that  spontaneous  magnetization  occurs  for  (\    >  ^ 


11a 


1.00 


*7K.  .50  - 


Fig. 7-  Isotherms  of  the  three-dimensional  spherical  model  of  a  ferromagnet 

based  on  equations  (38)  -  (i+O)  .  The  normalized  magnetization  M   =  M/|i 
is  shown  as  a  function  of  the  normalized  externa!!,  field  f/    =  H/kT  for 
various  constant  values  of  v'  =  j/kT.   Spontaneous  magnetization  occurs 
/^  =  .2527.   This  is  shown  by  the  fact  that  for  such  values  of 


for 


,  the  magnetization  is  different  from  zero  for  ?/  =  Q. 
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In  Appendix  D  it  is  proved  that 


dWI 
d?f 


=  1  (1  =  0  (42) 

^--0  ' 


^Uo^^F^  °^^^^^ 


d^ 

d?^ 


y=o 


(kk) 


Here  z  is  the  solution  of  equation  (k)   of  Appendix  C.  Magnetization  curves  based 
on  (38)  -  (4l)  are  shown  in  Figure  7.  The  spontaneous  magnetization  curve  based 
on  (^Ic)  is  shown  in  Figure  8. 

The  equation  of  state  of  the  three  dimensional  spherical  lattice  gas  is 
now  obtained  by  setting  n  =  J  in  (ll)  and  (15)  which  become 

Here  f,(z)  is  given  by 

2ji:  27t  2n 

f,(z)  =  i-  dw  dw  dw  In  (z  -  cos  w  -  cos  w  -  cos  w  )  .     (i+T) 

5      {2^y      j   i    ^    1   2  5  1        2        3 

^'000 


Equations  (12)  and  (59)  yield 


(W) 


12a 


0.0 


0      0.2     0.4     0.6     0.8      10 


^y, 


Fig. 8-  Spontaneous  magnetization  of  the  three-dimensional  spherical  model  of  a 
ferromagnet  based  on  eqtiation  (i+lc).  The  ordinate  is  the  limiting  value 
of  the  magnetization  Vl   (  /V,  (^  )  as  jV  tends  to  zero  through  positive  values, 
It  is  just  the  intercept  of  the  isotherms  with  the  axis  ^  =  0  in  Figure  ?• 
The  abcsissa  is  f   /  (j  - 


-  13 


From  (41),  (1+5)  and  (48)  we  derive,  in  Appendix  C,  the  entries  of  Table  k. 

In  this  table  z  is  the  solution  of  equation  (4)  of  Appendix  C,   while  v  and  v 
'■'  L      G 

are  determined  by  the  equations 
-1 


'  1  +  M(0+,  y ) 


^  .^  =  1  . 
^L    ^G 

For  very  small  and  very  large  v,  the  isotherms  behave  as  they  do  in  the  one- 
and  two-dimensional  cases.   However  for  0    >    ^  )    a   discontinuity  occurs  in 
the  p  ,  V  curves  at  the  parameter  value  V  =  0.  At  this  value  we  obtain  two 
different  specific  volumes  for  each  p  .   We  designate  the  smaller  v  to 
indicate  the  value  at  which  the  liquid  begins  to  vaporize.   The  larger  we  call 
V  ,  which  is  the  specific  volume  at  which  the  gas  begins  to  liquify.   The  jump 
in  specific  volume  is 

G    L   1.^2(0+,  rt) 


We  now  eliminate  y  between  (4-5)  and  (4-9),  obtaining 


*  2 

P   = 


^' 


2 


ml^j   .ln|l-  -f^(5) 


Equation  (52)  defines  that  half  of  the  phase  transition  ciirve,  or  boundary  of 

the  two  phase  region,  in  the  p  ,v  plane  at  which  the  liquid  begins  to  vaporize. 

The  other  half  of  this  curve  is  obtained  by  eliminating  v  from  (52)  by  means 

of  (50),  which  yields  a  relation  between  p  and  v_ .   This  is  the  half  of  the 

G 

* 
phase  transition  curve  at  which  the  gas  begins  to  liquify.   The  maximiun  of  p 

on  this  curve  occurs  at  v^.  =  v  =  2,  where  Sp  /5v^  =  3p  /Sv_  =  0.   It  is 
li    (j  li         G 


(^9) 


(50) 


(51) 


1  <  v^  <  2  .        (52) 


-  15a 


V  V 


o  + 


^h 


H^ 


^^ 
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interesting  to  note  that  at  the  maximum  5  p  /Sv_  =  S  p  /Sv.,  j-  0.  Thus  the 
phase  transition  curve  is  not  as  flat  at  its  maximum  as  is  that  of  the  two 
dimensional  ordinary  lattice  gas.  The  phase  transition  curve  or  two  phase 
boundary  is  shown  as  a  dashed  line  in  Figure  9- 

In  order  to  examine  the  isotherms  of  the  three  dimensional  spherical  lattice 
gas,  we  prove  in  Appendix  D  that 


d^ 
dv 


d^ 
dv 


^^2  =  -{V5)  /  <  U  <53) 


I    >    h-  (5^> 


From  (5^)  we  see  that  the  isotherms  meet  the  two-phase  boundary  with  zero 
slope  and  that  the  compressibility  is  infinite  at  this  boundary.  The  isotherms 
are  shown  in  Figure  9-   The  horizontal  line  within  the  two  phase  region  in  that 
figure  represents  the  discontinuity  (51)- 

Although  Figure  9  bears  a  marked  qualitative  resemblance  to  the  p,  v 
diagram  of  a  real  gas-liquid  system,  it  is  unrealistic  for  large  ^   at  small 
V  and  for  all  fl  at  large  v.   This  is  shown  by  the  small  negative  trough  of 
the  condensation  curve  for  1  <  v  <  2  and  by  a  corresponding  one  of  the  same 
depth  but  of  infinite  extent  for  large  values  of  v. 

6 .   Final  remarks . 


Our  method  of  computation  of  the  isotherms  from  (^5)  -  (48)  is  described 
in  Appendix  A.   In  this  computation  for  each  0'  ,  values  were  assigned  to  the 
two  parameters  z  and  "jf  ,    and  corresponding  values  of  p  and  v  were  computed. 
The  values  of  ^f"  ranged  from  -oo  to  +00  and  those  of  z  ranged  from  z  to  +oo  . 
The  minimiim  value  z   is  determined  in  Appendix  C  as  a  root  of  equation  (l)  of 
that  Appendix.   We  now  wish  to  determine  whether  the  isotherms  can  be  continued 
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beyond  the  two  phase  boundary  by  choosing  z  <  z  .   If  so^  we  wish  to  see 
if  such  a  continuation  has  physical  significance.   In  particular,  does  it 
yield  metastable  states? 

We  first  consider  the  spherical  Ising  model  and  write  (l6)  and  (l8)  in 
the  form 

w'  =  1  -  r^  ar  '  (55) 

V-  =  2^(J  -  n)n,  ■  (56) 

In  one,  two  and  three  dimensions  the  following  remarks  apply.   If  n  <  z  <  z 

1   -^^n^^  ^ 
it  follows  that  777  — >  1.   Then  (55)  shows  that  ?7J.  is  imaginary.   If 

^r     ^^  df^ 

|z|  <  n  then  from  (17)  we  find  that  - —  is  infinite.   Hence  we  need  consider 

only  z  <  -n. 

df^(z) 
From  (17)  we  note  that  — is  an  odd  f\inction  of  z.  Therefore  we  replace 

z  by  -z  in  ("55)  and  (56)  and  obtain 

y-    =   -2  ^m^   +  n),  z  >  n  .  (58) 

df 

Since  t^   is  finite  and  positive  for  z  >  n,  (57)  and  (58)  yield  two  real 
dz 

continuations  for  JH   and  ^^  since  two  values  of  the  square  root  can  be  used  in 
(57).  According  to  (58),^  and  ^  have  opposite  signs  for  each  solution,  and 
from  (57)  |7/l/|  >  !•   Thus  both  of  these  continuations  seem  to  be  physically 
unrealistic . 

We  shall  now  consider  the  possibility  of  continuing  the  isotherms  of  the 
spherical  lattice  gas.   First  we  consider  the  range  n  <  z  <  z  .   Since  "UJ/   is 
imaginary  in  this  range,  we  see  from  (12)  that  v  will  be  complex.   On  the  other 
hand,  if  z  <  n,  (17)  shows  that  f  (z)  is  complex.  Then  (I5)  and  (ll)  show  that 
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p   is  not  real.   Thus  if  z  and  )r  are  real,  there  is  no  real  continuation  of  the 

isotherms . 

Next  we  shall  show  that  when  v,  p*  and  (/  are  all  real  then  z  and  M    must  be 

real.   To  show  this  we  first  note  that  from  (l2).  Til    is  real  if  v  is  real.   Then 

from  (55)  it  follows  that  z  is  real.   For  if  z  were  not  real,  by  (I7)  df  /dz  would 

not  be  real  and  then  (55)  would  show  that  /?/  would  not  be  real.   Finally  from 

(18),  the  reality  of  z  and  ^' imply  that  ^^  is  real.   This  result  and  that  of  the 

preceding  paragraph  show  that,  on  the  basis  of  our  equations,  there  is  no  real 

continuation  of  the  isotherms  of  the  spherical  lattice  gas  into  the  two  phase 

region.   Thus  the  boundary  of  the  two  phase  region  is  the  natural  boixndary  of 

the  function  p*(v,  (1  )  . 

Appendix  A 

Method  of  computation. 

Each  isotherm  of  the  spherical  model  of  ferromagnetism  and  of  the  spherical 
lattice  gas  was  computed  by  assigning  a  value  to  Q   and  selecting  a  set  of  values 
for  z.  For  these  values  7II,')4,y,Y'   and  p*  were  computed  successively  from  (55), 
(56),  (12),  (15)  and  (11)  using  a  Burroughs  220  electronic  computer.   The  only 
difficulty  was  that  of  evaluating  f  (z)  and  df  (z)/dz  for  n  =  2  and  3.  The 
second  term  on  the  right  side  of  (31),  which  gives  df  /dz,  involves  the  complete 
elliptic  integral  of  the  first  kind.  This  Integral  also  occurs  in  (39)  which 
gives  df  /dz.   It  was  approximated  by  a  formula  of   C.  Hastings'-  '    ^'         -^  .     Then 
df  /dz  was  computed  by  numerical  integration  using  Simpson's  one-third  rule,  which 
was  also  used  to  compute  fp(z)  from  (56). 

The  fxinction  f,(z)  was  computed  from  the  series 
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5  kz'^  32z^     32z^     l,02l+z^ 


This  series  can  be  derived  by  first  writing  (^7)  in  the  for 


f^(z)  =  m  z  -  f:   ^  (A. 2) 


The  coefficients   in   (A. 2)   are  found  to  be 

«  n   It       ,  '  ^N^  n 

A^=^       III      '   E     cosco.   j        dc^dc^du,^    .  (A. 5) 

"      o  o  o    V"'  y 

The  integral  in  (A. 5)  is  expressible  in  terms  of  I  ,  the  modified  Bessel  fimction 
of  the  first  kind  of  order  zero.   This  yields 


a-*0   da     V       / 

Since  I  (a)  is  an  even  function  of  a,  it  follows  from  (AA)  that  A  =  0  when  n 
o^  '  ^    '       n 

is  odd.   Upon  substitution  of  the  power  series  for  I  into  (A.'^)  we  obtain  for 
even  n  =  2s 

A      (2s  -  1):   s  k 

'^  Z   [(s-k):]-2  ^  [o:(k-j).-]-2  (A.5) 

If  A   is  computed  for  s  =  1,  2,  5>  ^  from  (A.5)  and  the  resiilts  are  inserted 
into  (A. 2),  the  result  is  (A.l). 


-  18  - 

Appendix  B 

Limit  properties  of  the  three-dimensional  spherical  magnet 
and  the  spherical  lattice  gas. 

We  shall  now  prove  various  statements  made  in  the  text  about  the  behavior 

of  the  equations  of  state  of  the  three  dimensional  spherical  models  of  the 

magnet  and  the  lattice  gas.   Let  us  begin  with  (^1).   From  (58)  it  follows 

that  z  -»  00  when  ^  ^  0.   Then  from  (58)  -  (h-O)   we  find  that  as  (I  -^  0, 


z 

202. 

We  now  solve  (B.l)  for  (/z  and  substitute  the  result  into  (B.2),  which  yields 
(4la). 

To  prove  (^Ib)  we  first  note  from  (38)  that  z  ->■  3  when  fl  -»  00.  We  also 
note  that  df  (3)/dz  =  2  Q  ,  which  was  proved  by  G.  N.  WatsonL^J.  Then  for  Q 
large,  (38)  can  be  written  as 


(B.l) 


%-^   •  (B.2) 


y- 


?y(z-5) 


/  c 

Or 


Thus  11m     Jl         A      =  1  which  completes  the  proof.   To  prove  (4ld)  we  see 
from  (17)  that  df  /dz  decreases  monotonically  from  2  Q,     at  z  =  3  to  zero  at 
z  =  00  and  that  it  behaves  like  z   for  large  z .   If  ^  <  ^   and  }f  -*  0,    it 
follows  from  (38)  that  z  ->  z  >  3,  where  z  satisfies  (C'l+J,  and  (^Id)  follows. 
If    11    >    ^      then  z  -^  3  in  (38)  when  ^  -*  0  and  (38)  becomes  (B.3)^  proving 
(4lc).   To  prove  (^le)  we  observe  from  (38)  that  when  V  ^  ^  ,  z  -^  '=°     and 
consequently  df  /dz  ->■  0.   Then 


(B.3) 
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This  is  (^le)  and  completes  the  proof  of  (^l) . 

We  shall  now  derive  the  values  of  p*  in  Table  k .      When  ?/  ^  +  «>  it 
follows  that  M  -*  ±  1,    since  ?ll    is  odd.   Then  from  (B.J+),z  ^  \]4\/2Q-^   «> . 
For  large  z,  f,(z)  '^  In  z  =  In  \y\/2  0,  .      Substitution  of  these  limiting  values 
for  z  and  f  (z)  into  (46)  as  |  ?/  |  ^  oo  yields 


1     m 


Fj---  I  m^TT  +  \y\  •  (B-5) 


Substitution  of  (B.5)  into  (45)  yields,  when  ^j/  -*  +  «>,  the  entry  on  line  1, 
column  k   of  Table  k .     When  n   ->■  -  ^   this  substitution  yields  the  entry  on 
line  k,    column  k   of  Table  k.     When  y-  -*  +   0    and  (^<j?,z^z  >5  and 
E/2  0  {z-3)   -*  0.   Substitution  of  these  values  into  (^5)  yields  the  entry 
on  line  2,    colimin  k   of  Table  k.      Finally  when  ^  -*  +  0  and  (J    >    ^    , 
z  ^  z  =3,     y/2  0{z-3)  ^±7n/{0+,Q).     Again,  substituting  these  values 
into  (45)  yields  line  3  of  Table  h. 

Appendix  C 
Determination  of  z  . 


We  know  from  (17)  that  df  (z)/dz  is  a  positive  monotonic  decreasing  function 
of  z  behaving  like  l/z  for  large  z.  We  also  know  from  (24)  that  df,  (l)/dz  =  0°, 

from  (51)  that  df  (2)/dz  =  »  and  from  Appendix  B  that  df  (5)/dz  =  2  Q    .      Therefore 

2 
z  corresponds  to  the  smallest  non-negative  value  of  Jfl'   that  can  be  realized  from 

(55).   For  n  =  1,  2  and  all  0  and  for  n  =  3,  Q    <  ^c  ""^^^^  ^^  "^^   "  °'  ^^^^^ 


^J^J 


n  m 
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=      2(!     .  (C.l) 


More  specifically 


(1+^)  n  =  l  (C.2) 


:     =  K     K(  — )  n  =  2  (C.3) 


'^^3(^m) 


dz        /  -^^  /    /  c 


From  Appendix  B  we  also  have 


„  =  5,  c  <  r  .  (C.4) 


r.-3.f>l^.  (C.5) 


Appendix  D 
Shapes  of  the  isotherms  of  the  spherical  lattice  gas. 

The  slope  of  the  isotherms  of  the  spherical  lattice  gas  can  be  expressed 
simply  in  terms  of  the  properties  of  the  spherical  model  of  a  ferromagnet.  To 
show  this  we  differentiate  (ll)  with  respect  to  v  and  (12)  with  respect  to  ff  . 
Upon  combining  these  two  equations  we  obtain,  after  use  of  (l^). 


dE^  ^   .   (1+^)  (D.l) 

'^   "    2^   ^ 

This  relationship  also  relates  the  ordinary  lattice  gas  and  the  Ising  model. 

Let  us  now  evaluate  dW/d^  for  the  spherical  model.  To  do  so  we 
differentiate  (l8)  and  (l6)  with  respect  to  }i    and  combine  the  two  results. 
We  then  obtain  for  the  spherical  magnet 
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dlTl 

1                 y 

o 
dz" 

d'H 

-   2y  (z-n) 

h  iJTti^      d^f 

Since  — ;r^  <  0  it  follows  from  (D.2)  that  ^  >  0  and  from  (D.l)  that 
,  c    ~  (X  n     — 

dz 


dE^ 


<  0. 


<i%(z) 


From  this  equation  and  (17)  it  follows  that  p —  increases  monotonically 

dz 
from  -00  at  z  =  5  to  zero  at  z  =  «>. 


If  0  <  0  < 


rt   then,  as  is  shown  in  Appendix  C,  when  ?r  ->  0  it  follows 


that  z  -^  z^  >  5  and  ^  ^  0.   Hence  d^f  (z)/dz^  -*  d^f  (z^)/dz^  ^   0.   Upon 
substitution  into  (D.2)  of  the  limits  attained  by  z,  Iff}/   and  d  f,(z)/dz  as 
}i  -^  0  \re   obtain  {k^) .      If    0    >    f      then  from  Appendix  C  we  find  that  as 
'H    -*  0,    z  ->  3  and  7}l^  ^  (l  -  fjC   )  >  0.   Hence  from  (D.5),d^f  /dz^  -* 

-   2^/^rt  (z-j)"""'^    .   Substitution  of  these  limiting  values  of  z,  ^     and 

d^f  (z)/dz^  into  (D.2)  yields  (4U). 

Equation  (55)  then  follows  upon  substitution  of  (^Id)  and  (^5)  into  (D.l), 
and  (5^)  follows  upon  substitution  of  (4lc)  and  (kk)    into  (D.l). 


(D.2) 


dv 

Let  us  now  deduce  (ii2)  -  (kk) ,    (55)  and  (5^).   If  j?  =  0,  (kz)   follows 
directly  from  (J+la).   For  0    >  0  ve   have  for  z  '-^  5  [2,  p.  585] 

-  P  ,  (z-5)V2]-'  .  (B.5) 
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